hfAgeHfHf2R Documentation

Plotting Hf growth lines

Description

Plots two-stage Hf growth curves in the binary diagram age vs. initial 176Hf/177Hf.

Usage

   
hfAgeHfHf2(which = rownames(Hfinit), RCC = R.crust, 
    xmin = 0, xmax = NULL, ymin = NULL, ymax = DM[,"176Hf/177Hf"],  
    evol.lines1 = FALSE, evol.lines2 = TRUE, violins = FALSE,
    age.hist = FALSE, kde.bandwidth = 30, rugs = TRUE, 
    pch = labels[which, "Symbol"], col = labels[which,"Colour"], 
    cex = labels[which, "Size"], new = TRUE, 
    main = "Two-stage Hf isotopic growth diagram",...)

Arguments

which

character or numeric; names or sequence numbers of the samples to be used

RCC

numeric; the 176Lu/177Hf ratio of the intermediate crustal reservoir for calculation of the two-stage model

xmin

numeric; minimum of the age axis in Ga

xmax

numeric; maximum of the age axis in Ga

ymin

numeric; minimum of the 176Hf/177Hf axis

ymax

numeric; maximum of the 176Hf/177Hf axis

evol.lines1

logical; plot the evolution lines from present to the 'age'?

evol.lines2

logical; plot the evolution lines from the 'age' to the DM?

violins

logical; overplot the violin plots for each group?

age.hist

logical; overplot the frequency histogram of U–Pb ages?

kde.bandwidth

numeric; KDE bandwidth.

rugs

logical; show the rugs on both axes?

pch

plotting symbols

col

plotting colours

cex

relative size of the plotting symbols

new

logical; should be a new plotting window opened?

main

main title for the plot

...

optional parameters to the underlying function plotWithLimits

Details

In this case, Hf growth curves in individual samples are plotted solely using a two-stage model (Kemp & Hawkesworth 2003, 2014). Moreover, shown are Hf growth curves for the two main mantle reservoirs, CHUR and Depleted Mantle (DM). The optional small ticks, or rug, on the x axis correspond to two-stage Hf model ages. On the y axis, the extra ticks portray the initial 176Hf/177Hf ratios.

Another option, controlled by the argument age.hist, is plotting of the kernel density estimate (KDE) curve for the U–Pb ages along the x axis. The smoothness of this curve can be controlled by the parameter kde.bandwidth, with the 'convenient compromise' σ = 30 (Vermeesch & Garzanti 2015; Andersen et al. 2018) set as a default value. In general, the higher is this value, the smoother curve results.

hfAgeHfHf1.png

hfAgeHfHf2.png

Value

None.

Plugin

Hf.r

Note

This function is Figaro compatible (except violins). For the Hf growth diagram of age vs. ε(Hf) values see hfAgeEps.

Author(s)

Vojtěch Janoušek, vojtech.janousek@geology.cz

References

Andersen T, Kristoffersen M, Elburg MA (2018) Visualizing, interpreting and comparing detrital zircon age and Hf isotope data in basin analysis – a graphical approach. Basin Res 30:132–147 doi: 10.1111/bre.12245

Blichert-Toft J, Albaréde F (1997) The Lu–Hf isotope geochemistry of chondrites and the evolution of the mantle–crust system. Earth Planet Sci Lett 148:243-258. doi: 10.1016/S0012-821X(97)00040-X

Bouvier A, Vervoort JD, Patchett PJ (2008) The Lu–Hf and Sm–Nd isotopic composition of CHUR: constraints from unequilibrated chondrites and implications for the bulk composition of terrestrial planets. Earth Planet Sci Lett 273:48-57. doi: 10.1016/j.epsl.2008.06.010

Chauvel C, Garçon M, Bureau S, Besnault A, Jahn BM, Ding ZL (2014) Constraints from loess on the Hf–Nd isotopic composition of the upper continental crust. Earth Planet Sci Lett 388:48-58. doi: 10.1016/j.epsl.2013.11.045

Griffin WL, Wang X, Jackson SE, Pearson NJ, O'Reilly SY, Xu X, Zhou X (2002) Zircon chemistry and magma mixing, SE China: in-situ analysis of Hf isotopes, Tonglu and Pingtan igneous complexes. Lithos 61:237-269. doi: 10.1016/S0024-4937(02)00082-8

Kemp AIS, Hawkesworth CJ (2003) Granitic perspectives on the generation and secular evolution of the continental crust. In: Holland HD, Turekian KK (eds) Treatise on Geochemistry vol. 3, The Crust (ed. R.L. Rudnick). Elsevier-Pergamon, Oxford, pp 349-410. doi: 10.1016/B0-08-043751-6/03027-9

Kemp AIS, Hawkesworth CJ (2014) Growth and differentiation of the continental crust from isotope studies of accessory minerals. In: Holland HD, Turekian KK (eds) Treatise on Geochemistry (Second Edition). Elsevier, Oxford, pp 379-421. doi: 10.1016/B978-0-08-095975-7.00312-0

Lancaster PJ, Storey CD, Hawkesworth CJ, Dhuime B (2011) Understanding the roles of crustal growth and preservation in the detrital zircon record. Earth Planet Sci Lett 305:405-412. doi: 10.1016/j.epsl.2011.03.022

Spencer CJ, Kirkland CL, Roberts NMW, Evans NJ, Liebmann J (2020) Strategies towards robust interpretations of in situ zircon Lu–Hf isotope analyses. Geosci Front 11:843-853. doi: 10.1016/j.gsf.2019.09.004

Vermeesch P, Garzanti E (2015) Making geological sense of 'Big Data' in sedimentary provenance analysis. Chem Geol 409:20–27 doi: 10.1016/j.chemgeo.2015.05.004

See Also

hfAgeEps1, hfAgeEps2, hfIso, hfDMage, hfDM2stgAge, hfViolinplot, vioplot

The actual plotting is done by the function plotWithLimits.

Examples

sampleDataset("khantaishir_Hf")
    
hfAgeHfHf2(1:39, RCC = 0.022, # Basaltic (Lancaster et al. 2011)
    evol.lines1 = FALSE, evol.lines2 = TRUE, xmax = 1.5, ymin = 0.2822, col = "blue" , cex = 1.5)

hfAgeHfHf2(1:39, cex = 1.5, xmax = 1.3, ymin = 0.2822, evol.lines1 = TRUE, evol.lines2 = FALSE) 

hfAgeHfHf2(1:39, xmin = 0.35, xmax = 1.3, ymin = 0.2822, ymax = 0.2830,
    evol.lines1 = FALSE, evol.lines2 = TRUE, age.hist = TRUE,
    kde.bandwidth = 48, # 'optimal' sigma value (Andersen et al. 2018)
    cex = 1.5, rugs = TRUE) 
        
hfAgeHfHf2(1:39, xmin = 0.35, xmax = 1.3, ymin = 0.2822, cex = 1, 
    evol.lines1 = FALSE, evol.lines2 = TRUE, violins = TRUE)